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ABSTRACT 



Context. The optically invisible infrared-source NGC 2264 IRS 1 lying north of the Cone Nebula is thought to be a massive young 
stellar object (~ 10 Mq). Although strong infrared excess clearly shows that the central object is surrounded by large amounts of 
circumstellar material, no information about the spatial distribution of this circumstellar material has been available until now. 
Aims. We used the ESO Very Large Telescope Interferometer to perform long-baseline interferometric observations of 
NGC 2264 IRS 1 in the mid-infrared regime. Our observations resolve the circumstellar material around NGC 2264 IRS 1, provide 
the first direct measurement of the angular size of the mid-infrared emission, and yield direct constraints on the spatial distribution of 
the dust. 

Methods. We analyze the spectrally dispersed interferometric data taken with MIDI at two different position angles and baseline 
lengths. We use different approaches (a geometrical model, a temperature-gradient model, and radiative transfer models) to jointly 
model the observed interferometric visibilities and the spectral energy distribution. 

Results. The derived visibility values between ~ 0.02 and ~ 0.3 show that the mid-infrared emission is clearly resolved. The char- 
acteristic size of the MIR-emission region is ~ 30 - 60 AU; this value is typcial for other YSOs with similar or somewhat lower 
luminosities. A comparison of the sizes for the two position angles shows a significant elongation of the dust distribution. Simple 
spherical envelope models are therefore inconistent with the data. The radiative transfer modeling of our data suggests that we ob- 
serve a geometrically thin and optically thick circumstellar disk with a mass of about 0. 1 Mq. 

Conclusions. Our modeling indicates that NGC 2264 IRS 1 is surrounded by a flat circumstellar disk that has properties similar to 
disks typically found around lower-mass young stellar objects. This result supports the assumption that massive young stellar objects 
form via accretion from circumstellar disks. 

Key words, techniques: interferometric - stars: individual: NGC 2264 IRS 1 - stars: formation-stars: circumstellar matter 



1. Introduction 

The question of how massive stars (> ^ Mq) form still re- 
mains unclear because observations of massive young stellar ob- 
jects (MYSOs) are rare (Zinnecker & Yorke 2007). This has sev- 
eral reasons. First of all, these objects are surrounded by large 
amounts of gas and dust, which absorb nearly all the optical and 
near-infrared (NIR) radiation from the central source, thus mak- 
ing a direct observation very difficult. Second, their large dis- 
tances mean that to study them we need data of high angular 

' resolution. Other difficulties are the fast evolution of MYSOs 

' and that they are rarely found in isolation. 

Two different formation scenarios for MYSOs are cur- 
rently discussed: the coalescence scenario and a formation sim- 
ilar to low mass protostars. In the first scenario, several low 
mass stars merge and finally form a single, rn ore massive 
star (lBonnelletalJll998l: IZinnecker & Batell2002h . In the sec- 
ond scenario, the radiation pressure problem that would pre- 
vent spherical accretion as soon as the protostellar mass ex- 
ceeds ~ IOMq is solved by stellar outflows, which open a 
cavity in polar direction, through which the radiation can es- 



* Based on observations collected at the European Organisation for 
Astronomical Research in the Southern Hemisphere, Chile, observing 
programs 076.C-0725(B) and 082.C-899(A). 



cape and the radiation pressure in the equatorial plane is low- 
ered. Therefore, an infall of material in this plane is possible. 
Observations have provided some evidence that massive stars 
can form in an analogous way to low mass star s. For instance, 
circumsteUar disks (e.g., IRAS 20126-1- 4104, ICesaroni etM 
200 51; AFGL 49 0, ISchrever et all l2006t IRAS 13481-6124, 
Kraus et al. '2010"), molecular outflows, or collimated jets (e.g., 
Cunningha m et al. 2009) are often interpreted as evidence of 
ongoing accreti on. From the theoretica l side, both the compet- 
itive a ccretion ('Bonnell & Bate! l2006l) and the turbulent core 
model dMcKe e & Tan 2003) agree with the formation of disks. 
Furthermore, simulations demonstrate that high-mass stars may 
form by disk accretion (Kuiper et al. 2010) . 

Traditionally, the spatial distribution of the surrounding dust 
was usually studied by the modeling the spectral energy distri- 
bution (hereafter SED). Such SED model fits are highly am- 
biguous (Thammetal. 1994; Men'shchikov & Henning 199"^. 
Therefore, more information is required to determine the real 
geometry of the circumstellar material. To test the existence 
of a circumstellar disk, and thus the disk accretion scenario of 
massive star formation, information at high angular resolution 
(coiTesponding to physical scales below 100 AU) is needed. 
Furthermore, only radiation at long wavelengths (e.g., infrared, 
sub-mm) is able to escape from the dense dusty surroundings of 
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the star. Infrared long-baseline interferometry can thus provide 
direct spatial information on the required small scales. 

NGC 2264 is a young, star-forming cluster with an age of 
about 3 Myr located in the Mon OB 1 molecular cloud com- 
plex. Literature values for the distance of the cluster vary be- 
tween 400 pc and 100 pc, although here we adopt the most 
recent value of 913 pc (iBaxteret all 12009!) . NGC 2264 IRS 1, 
the brightest IR-source in this region, was discovered in 1972 
by D. Allen (hence is also called A llen's Source) and has 
no optical counterpart. lAUenl (Il972h pointed out that since 
NGC 2264 IRS 1 is the brightest and most luminous source in 
this field, it is the likely source of the radiation pressure creating 
the Cone Nebula. It is believ ed to be a young star wi th a mass of 
9.5 Mo, spectral ty pe B0-B2 dThompson et al.ll 19981) . a luminos- 



ity of 3.5 ■ 10^ Lq dHarvev et alj|1977 ). and a visual extinction 



of 20-30 mag. However, owing to its high extinction in the opti- 
cal, all values are rather uncertain and should be in terpreted with 
care, e.g., luminosity values between 2.3 ■ 1 0"^ L© (iNakano et al.l 
l2003h and 4.7 ■ 10^ Lq (ISchwartz et alJll985h can be found in the 
literature. 

NGC 2264 IRS 1 is associated with molecular outflows 
and dense molecular clumps, and in near-infrared images a 
twisted jet-like feature can be seen extending n orth of the source 
dSchreveretaLl 119971; IWard-Thompson et all [2000). The mor- 
phology of the circumstellar environment is unclear. Several 
studies did n o t find any hint of an asymmetric structure. 
ISchrever et al.l (l2003h were unable to detect any si gns of a cir- 
cumste llar disk in their molecular line study and Ide Wit et al.l 
(|2009|) do not discern any disk-like structure in their images 
taken in the MIR with COMICS. However, models of a spher- 
ical symmetric distribution of the circumstellar material fail 
to reproduce the SED of the o bject in its entirety and were 
ruled out in severa l publications ([Thompson & Tokunagal[l978t 
iHarvev et al.ll 19771) . However, the spatial resolution of all these 
studies may not be suflicient to resolve a disk-like structure. 

In this paper, we present mid-infrared interferometric obser- 
vations of the MYSO NGC 2264 IRS 1 performed with the Mid- 
Infrared Interferometric Instrument (MIDI) at the Very Large 
Telescope Interferometer (VLTI). The MIDI data and their re- 
duction are shown in Sect. |2] Diff'erent modeling approaches 
(geometrical, temperature-gradient, and radiative transfer mod- 
els) for the visibilities and the SED are presented in Sect. |3] 
Implications for the object and final conclusions are discussed 
in Sect. H 



2. Observations and data reduction 

NGC 2264 IRS 1 was observed on 2005 December 23 in the 
course of ESO observing programme 076.C-0725(B) (PI: S. 
Kraus) and on 2009 February 15 in the course of ESO pro- 
gramme 082.C-0899(A) (PI: M. Feldt) with MIDI at the VLTI 
(see Table |2]i. MIDI (Leinert et al. 2003) is a two-telescope 
beam-combiner working in the wavelength range from 8 to 
I3fim and using the principle of a Michelson interferometer. 
To measure the fringe amplitude, the optical path difference 
(OPD) is varied with internal delay lines. The 2005 data of 
NGC 2264 IRS 1 were taken with a baseline of 89.1 m in the 
SCI-PHOT mode, in which interferometric and photometric data 
are measured simultaneously. The grism (a combination of a 
grating and a prism) with a spectral resolution of ^ = 230 was 
used. In the same night, two calibrators were observed with the 
same configuration. HD 39425 a few hours earlier and HD 61421 
directly after the object. The 2009 data were taken with a base- 
line of 40.2 m in the HIGH-SENS mode, in which interferomet- 



ric and photometric data are obtained separately. For this obser- 
vation, the prism with a spectral resolution of ^ = 30 was used. 
HD 48217 and HD 49161 served as calibrators. 

For data reductio n and calibration, w e used the software 
package MIA-i-EWS dKohler & JaffdllOOSl) . With this software, 
the correlated flux, visibility, and the total flux can be deter- 
mined. Typical errors in one s ingle MIDI visibility measurement 
are about 10% dLeinert et al.l 200 4). As an estimate of the error 
bars, the standard deviation in the results for different calibrators 
was used. As the conditions during the observation in 2009 were 
not photometric, we did not use the 2009 photometry data, but 
the photometry from the 2005 measurement to obtain the vis- 
ibility. Figure [1] (left) shows the calibrated visibilities for both 
observations. 

Furthermore, a low-resolution IR spectrum of 
NGC 2264 IRS 1, measured with the SWS instrument 
aboard the ISO satelHte is available (TDT No. 71602619, PI: D. 
Whittet), as well as fluxes in several other wavelength bands. 
These data points are listed in Table [T| Fig. [T] (right) shows the 
total flux spectrum of NGC 2264 IRS 1 as measured by MIDI 
in 2005 compared to the total flux spectrum measured by the 
ISO satellite. Both spectra display a deep absorption feature 
over the whole MIDI wavelength range caused by silicates. 
The shape of both spectra (MIDI and ISO) is similar but the 
flux level of the MIDI spectrum is about 25 Jy (about 20 %) 
lower, which might be due to the different beam sizes of the 
instruments (see Table [1). The spurious feature in the MIDI 
spectrum at wavelengths between 9.3 yum and 9.8 yum is due to 
atmospheric ozone. Se veral absorption features can be identified 
in the ISO spectrum (Gibb et al] 12004 ). The broad absorptions 
around 3.1 and 6.0 /im arise mainly from water ice, while the 
carrier of the 6.8 //m feature is more uncertain, but might be 
associated with solid NH|. A rather sharp feature at 4.27 //m 
is commonly attributed to CO2 ice (Guertleret al. 1996), as is 
the absorption seen at 15.2 yum. The feature around 4.7 jum is 
probably a combination of CO gas line absorption and a variety 
of ice absorption features. All these findings point to relatively 
large column densities of cold dense material along the line of 
sight towards the central source. 



3. IVIodeling 

We now attempt to derive the properties of NGC 2264 IRS 1 
from our measurements using different modeling approaches. 
First, we use simple geometrical models to characterize the 
size of the dust emission region (see Sect. 13. 11 1. To test for the 
possibility of a circumstellar disk, we compare the data with 
a temperature-gradient model for such a disk (Sect. I3.2l i. To 
check the parameter space over a wider range, w e employ the 
online model SED fitter by Robitaille et al.1 iJOQlh in Sect. [33] 



A more detailed modeling is done in Sect. 13.41 where the two- 
dimensional radiative transfer code RADMC is used. 



3. 1 . Simple geometrical models 

As a first step, one can derive an estimate of the diameter of the 
mid-infrared (MIR) emission region of the object from the visi- 
bilities by usin g simple geometrical model s. F or more details o f 
this we refer to lBerger & SegransanI d2007l) and lTristranil d2007l) . 
We first consider the possibility of a binary. Its existence 
would result in a sinusoidal variation in the visibility (ver- 
sus wavelength), ranging between a maximum and a minimum 
value. The maximum value depends on the extension of the 
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Table 1. Measured fluxes for NGC 2264 IRS 1 



Wavelength 


Flux 


Instrument/Reference 


Aperture 


[jum] 


[Jy] 






2-45 


9-670 


LSOfSloanetal.f20031> 


14" X 20" 


1.65 


0.9 ±0.04 


2-MASS 


4" 


2.2 


6.8 ±0.016 


2-MASS 


4" 


8-13 


60-130 


MIDI 


0.5" 


12.0 


146 ±9 


IRAS 


0.75' X 4.5' 


24.5 


330 


COMICS 


42" X 32" 


25.0 


324 ±13 


IRAS 


0.75' X 4.6' 


53.0 


980 ±50 


Harvev et al. (1977) 


17" 


60.0 


911 


IRAS 


1.5' X 4.7' 


70.0 


960 ±290 


Sargent etal. (1984) 


3' 


100.0 


1560 


IRAS 


3'x5' 


100.0 


1645 ±82 


Harvev et al. (1977) 


28" 


175.0 


1530 ±77 


Harvev etal. (1977) 


46" 


350.0 


188 ±70 


Chinietal. (1986) 


30" 


1300 


13 ±2.5 


Chini etal. (1986) 


90" 



Table 2. Observation log of the MIDI measurements of NGC 2264 IRS 1 



Date 



Proj. Baseline 
[m] 



PA 

n 



Mode 



12/23/2005 
02/15/2009 



89.1 (UT2-UT4) 

40.2 (UT2-UT3) 



81.1 

43.9 



SCI-PHOT&GRISM 
HIGH-SENS&PRISM 




1000 



100 



10 



X [urn] 




100 



Fig. 1. Left: Visibilities of NGC 2264 IRS 1 measured with MIDI. The visibility points measured in 2005 (projected baseline 
length 89.1 m) are shown in red, those measured in 2009 (projected baseline 40.2 m) are shown in green. Right: Spectrum of 
NGC 2264 IRS 1 as measured with MIDI (red crosses) in 2005 compared to the ISO spectrum (blue line). 



MIR emission regions (and would be 1 for point like sources), 
whereas the amplitude depends on the flux ratio of both compo- 
nents (the minimum would be for a flux ratio of 1:1). Owing 
to the slit width of MIDI (see Table [T]), the maximum separa- 
tion of a binary that can be identified with this instrument is 
~ 230 AU. The minimum separation at which a binary could 
be resolved with the 90 m baseline is ~ 30 AU (diffraction limit; 
half a period of the modulation in the visibility would be seen). 
Both observed visibilities show no typical sinusoidal variation, 
hence a binary system with a separation of > 30 AU is unlikely 
and not considered during the modeling process. However, with 
our measurements we could miss very faint binary companions 
or systems with separations of about 30 AU oriented non-parallel 
to the position angle of the 90 m baseline. 

We consider two different geometrical models for a single 
source: a Gaussian distribution and a uniform disk (UD). The 



results can be seen in Fig.|2] All curves show an increasing size 
of the object with increasing wavelength, as expected because of 
the different dust temperatures in different parts of the dust distri- 
bution. For the measurement performed with the 90 m baseline, 
the size varies between 16 and 30 mas for the Gaussian model 
and between 20 and 35 mas for the uniform disk model, whereas 
for the 40 m baseline, the size of the object lies between 30 and 
44 mas for the Gaussian model and between 44 and 64 mas for 
the uniform disk model. As the observations were taken at dif- 
ferent position angles, the diff'erent size estimations suggest an 
elongation in the dust distribution (e.g., a circumstellar disk or 
a flattened envelope). Although the intensity profile of the real 
source possibly differs from those distributions, the Gaussian 
model in particular is usually a good approximation of the mid- 
infrared morphology of YSOs. Despite all the uncertainties in 
the visibility measurements, hence the size estimation, the fitting 
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Fig. 2. Mid-infrared diameter of NGC 2264 IRS 1 as a function 
of wavelength for the 2005 data (crosses) measured with a 90 m 
baseline and the 2009 data (rectangulars) measured with a 40 m 
baseline. The values for the diameter in AU (right y-axis) were 
computed using a distance of 913 pc. Two different models were 
used: a Gaussian model (red) and a uniform disk model (green). 



of the Gaussian model still helps us identify an asymmetry with 
a significance of 3 cr and the UD model within 2 cr. This result 
clearly suggests that the MIR emission region is not spherically 
symmetric but elongated. 

3.2. Temperature gradient model 

Many young stars of both lo w and high mass are known to 
have circumstellar disks (e.g., iKrauset al]|2010t iLeinert et al.l 
l2004t iRatzka et al.ll2009l) . As a first approach to interpreting si- 
multaneously SED and N-band spatial data, we therefore use a 
te mperature-gradient m odel for a circumstellar disk as described 
in lMalbet et al.l (l2005h . The disk specified by this model is opti- 
cally thick, but geometrical thin with an inner radius r^in and an 
outer radius rmax- The distribution of the temperature in the disk 
is described following the equation 



(1) 



T(r) = To - 



where Tq is the temperature at a reference radius ro = 1 AU. 
The exponent q is expected to lie approximately between 0.5 for 
passive, i rradiated, flared dis ks and 0.75 for passive, irradiated, 
flat disks dMalbet et al.ll2005h . 

One can imagine such a disk to be composed of rings with 
radius r and thickness dr. If every ring emits blackbody radia- 
tion, one can evaluate the flux by integrating over the radius r 
such that 



^,i(0) 



2n r- 



rB;i(T(r)) cos(i)dr. 



(2) 



where d is the distance to the object, B^ the Planck function for 
a blackbody with temperature T(r), and i the inclination of the 
disk. For the visibility, one then gets 



b - Juq + Vq ■ cos^(i). 



Here, © is the position angle of the disk and m© and vq are given 
by 

uq - u cos - v sin 

and 

v© = M sin -H V cos 0. 

We calculated a large number (^ 1000) of temperature-gradient 
models and tried to find a model that is simultaneously able to re- 
produce all observations (both visibilities and SED). For this, we 
varied the innerradius between 0.5 and 2 AU (in steps of 0.5 AU) 
and Tq between 1400 and 3000 K in steps of 200 K. For the pa- 
rameter q, we tested the values 0.5, 0.6, and 0.7. Inclination and 
position angle were varied in steps of 10°. 

The temperature gradient model with the parameters Usted 
below is able to reproduce the SED in the mid- and far-infrared, 
but fails in the near-infrared. All of the calculated models failed 
completely to reproduce the shape of the visibility measured 
with the 40 m baseline (rising values going to larger wave- 
lengths). As the best-fit model, we therefore show the model 
that most successfully reproduces the other visibility curve and 
the SED simultaneously, although different models led to a bet- 
ter fit for the 40 m visibility at least with respect to its absolute 
level. The best-fit model is shown in Fig. [3] and has the following 
parameters: 

- r,ni„ = 0.5 AU 

- r^ax = 4000AU 

- To = 2350 K 

- i = 75° 

- ^ = 0.5 

- = 90° 

The spectral features seen in the SED and the effects of the sili- 
cate absorption on the visibility cannot be reproduced, as radia- 
tive transfer effects (emission, absorption, and scattering) are not 
included in the model. 

A temperature of To = 2350 K at a reference radius of 1 AU 
would imply a temperature of 3250 K at the inner disk radius 
of 0.5 AU for this model. This is much too hot for the existence 
of normal dust grains. A po ssible explanation could be the ex- 
istence of an inner gas disk dKraus et alj|2008h or of dust grains 
with higher sublimation temperature (iBenistv et alJl2010l) . 



3.3. Disk and envelope models in the Robitaille grid 

To check a wider parameter space for different models (espe- 
cially for the possibility of a circumstellar envelo pe), we em- 
ployed the online SED fitting tool developed by .Robitaille et al.l 
(l2007h . This tool compares 200,000 models of flie SEDs of 
YSOs, which were precomputed using a 2D radiative transfer 
code (Whitney et al. 2003), with the measurements. These mod- 
els are described by a set of 28 parameters, such as the proper- 
ties of the central object (e.g., mass, temperature), the proper- 
ties of the circumstellar envelope (e.g., outer radius, accretion 
rate, opening angle), the circumstellar disk (e.g., mass, outer 
disk radius, flaring), and the ambient density (for details see 
iRobitaille et al.ll2006h . In addition to defining SED points, the 
user has to set a range for the interstellar extinction and for the 
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Fig. 3. Left: Comparison of the best-fit temperature-gradient model to the visibilities of NGC 2264 IRS 1 with the model parameters 
described in Sect. l3.2l (green for the 40 m baseline, red for the 90 m baseline). The resulting modeled visibility for the 40 m baseline 
is shown by the gray line, the one for the 90 m baseline by the black line. Right: Comparison of the SED of our best-fit temperature- 
gradient model (black) described in Sect. 13. 21 with the ISO spectrum (blue) of NGC 2264 IRS 1. The MIDI spectrum is shown with 
red crosses, all other measurements from Table[T]are shown with black rectangulars. 



distance of the object. We chose the interstellar extinction Ay 
to be between and 40 mag and the distance d to be between 
800 pc and 1000 pc. 

The number of points that can be fitted by the SED online 
tool is limited. Therefore, from the ISO spectrum about one flux 
point per one micron was used. From the MIDI spectrum, about 
two flux points per micron were chosen as the shape of the spec- 
trum changes much more in this wavelength regime because of 
the silicate absorption feature. 

When comparing the model visibilities and fluxes to the ob- 
served data, it is important to consider the very different (effec- 
tive) beam-sizes of the different observations. The MIDI instru- 
ment has an effective field-of-view of 0.5" (slit width), corre- 
sponding to a! 450 AU at the distance of our target. The beam 
sizes for all other data points beyond near-infrared wavelengths 
are much larger: the beam size from which the ISO spectrum 
was extracted (* 13 000 x 18 000 AU) is already * 1 100 times 
larger, and the beam sizes for all far-infrared data points are at 
least ~ 3000 times larger than the MIDI field-of-view. The ISO 
and far-infrared fluxes do therefore not only trace the emission 
from the central YSO and its immediate circumstellar material 
(on spatial scales of up to ~ 500 AU), but also contain contribu- 
tions from the surrounding extended molecular cloud, on spatial 
scales of > 20000 AU (^> 0.1 pc). 

Therefore, as a first approach we include all SED points with 
A > 13//m only as upper limits. T o compute the visibilities, we 
used the HO-CHUNK code from I Whitney et al J ( I2b03i) to cal- 
culate images for the model. Doing so, the ten best-fit models 
consist of circumstellar disks only (see Sect. l3.3.TI ). As a second 
approach, none of the data points was used as an upper limit, 
but they all were given with their errors. Here, rather different 
best-fit models are found by the fitting tool. They all consist of a 
large circumstellar envelope and either have no additional disk, 
or just a small ring-like structure with radii from 0.5 to 8 AU (see 
Sect. l3.3.2l i. For similar appHcations of this fitting tool and its ad- 
vantages and lim i tations in this respec t, we refer to for instanc e 
iLinz et al.1 (l2009h . lde Wit et al.1 (l2010h . and lFoUert et aP (l2010h . 
We note that the parameter space is not covered uniformly. In 
particular, the spacing of the grid parameters is much finer in the 
case of low-mass YSOs than for MYSOs. 



3.3.1. Robitaille disk models 

The SED and visibilities of the best-fit Robitaille disk models 
are shown in Fig.|4] This model is able to reproduce the flux in 
the near- and mid-infrared. To reproduce the flux in the FIR, one 
would have to add for example some blackbody components as 
shown later for the RADMC model (see Sect l3.4l ). The model 
can reproduce the level of the 90 m baseline visibility, and also 
the 40 m baseUne visibility is reproduced within the errorbars. 
The parameters of this best-fit model are (model ID 3004478): 

- Mstai-= 11.6 Mo 

- r,tar = 27,720 K 

- Lstai-=1.21X1O4L0 

- rout = 84 AU 

- /idi.sk(100AU) = 5.6AU 

- mdisk = 0.09 Mq 

- / = 70° 

- Av = 20.27 mag 

- = 40° 

Here, rout is the outer disk radius and /idisk would be the height 
of the disk when extending it to a radius of 100 AU. 

3.3.2. Envelope plus disk models 

The Robitaille envelope model with the smallest x^ is shown in 
Fig- El It fails to model especially the shape of the 40 m baseline 
visibility curve. The flux in the mid-IR can be reproduced, but 
the model cannot reproduce the flux in the near-IR. To attain the 
flux level in the far-IR, one would once again have to add, e.g., a 
blackbody component. The parameters of the best-fit model are 
(model ID 3019090): 

- M3tat- = 8.5M0 

- r„ar = 4200K 

- ^star - 82 Rq 

- Lstar = 2.O5xlO3L0 

- /-envelope = 36,000 AU 

- '^envelope - 1.84 Mq 

- '•disk = 3.6 AU 

- /!disk(100AU) = 7.7AU 
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Fig. 4. Left: Comparison of the best-fit Robitaille disk model (black crosses 40 m baseline, black rectangulars 90 m baseline) for the 
visibilities of NGC 2264 IRS 1 with the MIDI visibilities (red 90 m baseline, green 40 m baseline). Right: Comparison of the best-fit 
Robitaille model (black) for the spectrum of NGC 2264 IRS 1 with the ISO spectrum (blue), the MIDI spectrum (red crosses), and 
the data points from Table [T] (black rectangulars). 



- mdisk = 0.2 Mq 

- i = 18° 

- Ay = 3.8 mag. 

A change in the position angle does not in this case lead to 
significant differences for this model, as the dust distribution is 
nearly symmetric. Therefore, a best-fit position angle cannot be 
given. As mentioned before, this model consists mainly of a cir- 
cumstellar envelope, with an outer radius of renveiope- However, 
there is also a disk component with a very small outer radius 
'"disk = 3.6 AU. The central object is much cooler and larger than 
the central object in the disk models. 

3.4. Two-dimensional radiative transfer modeling 

A more detailed physical model for a circumstellar disk is 
obtained by using the two-dimensiona l (2-D) radiative trans- 
fer m odel code RADMC (Version 3.1. iDullemond & Dominila 
12004 ). which provides calculations of continuum radiative trans- 
fer in three-dimensional axisymmetric (therefore 2-D) circum- 
stellar dust distributions around a central illuminating star. The 
main radiative transfer calculations are done with a Monte Carlo 
algorithm. To produce spectra and images, a post-processing 
with a ray-tracing code is done. The disk is described us- 
ing a passive irradiated fl ared disk based on the model of 
IChiang & Goldreich'(1997). but the code also allows the user to 
insert a vertical pufFed-up inner rim to produce a self-shadowed 
disk. The space is mapped in spherical coordinates. We therefore 
have radial grid points r and angular grid points 0, where 6-0 
is the polar axis and = flmax = f is the equatorial plane. We 
again tried to reproduce the SED and visibilities simultaneously 
and also to fit the silicate absorption feature. To scan the param- 
eter space, more than 10,000 different models were calculated 
and compared to the observed data. 

The geometry of the disk is described by the pressure scale 
height /idisk at a reference radius rjisk, and an exponent p in the 
following way (for ri,, < r < rdisk): 

{-)'■ 

\'"disk/ 



h(r) = h 



disk 



(4) 



Furthermore, either the surface density at the outer radius ctq can 
be given by the user directly, or (as we did) be calculated from 



the user-defined total disk mass mdisk = ^Wgas + '"dust- We used 
a gas-to-dust ratio of 100. The surface density as a function of 
radius is then given by 



cr(r) = CTQ ■ 

V '"disk 



(5) 



For r > rdisk, the exponent q is chosen to be equal to 12 such 
that the density drops so fast that rdisk is eff'ectively the outer 
disk radius. This outer radius was varied between 30 AU and 
150 AU. The density distribution in the whole disk is described 

by 



p(r,0)' 



h(r) ■ r 



■ exp 



1 



h(r) 



(6) 



Parameters that remained fixed during the modeling process, 
were the temperature of the central object Tstar - 25, 000 K (ac- 
cording to the spectral type B0-B2), the exponent ^ = 12 for 
f > /"disk, and the temperature at the inner rim Ti„ = 1700 K. 
For the chemical compositio n of the dust, the o ptical proper- 
ties of astronomical silicate (iDraine & Leelll984h and graphite 
were chosen, u sing the grain size distribution of the MRN model 
(lMathisetal.11977.) . A mixture of 70% sihcate and 30% graphite 
was assumed. 

In addition to the original RADMC code, we introduced 
foreground extinction Ay, as otherwise all models failed to re- 
produce the deep silicate absorption feature. Parameters that 
were varied during the fitting process are shown in the following 
list together with their values for the best-fit RADMC model: 

- '"disk = 50 AU 

- /^disk = 0. 1 

- Lstar = 3.98x 10^ L© 

- p = 2/7 

- '"disk = 0.1 Mo 

- ^ = 2.5 for r < rdisk 

- / = 30° 

- PA = 40° 

- Ay = 27 mag. 

This best-fit model is shown in Fig. |6] It provides a good fit to 
the visibilities as well as the SED up to 12 /im. 
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Fig. 5. Left: Comparison of the best-fit Robitaille envelope model ( black crosses 40 m baseline, black rectangulars 90 m baseline) 
for the visibilities of NGC 2264 IRS 1 with the MIDI visibilities (red 90 m baseline, green 40 m baseline). Right: Comparison of 
the best-fit Robitaille model (black) for the spectrum of NGC 2264 IRS 1 with the ISO specti'um (blue), the MIDI spectrum (red 
crosses), and the data points from Table[T](black rectangulars). 



The RADMC model describes the central YSO and its cir- 
cumstellar matter up to the outer disk radius of 50 AU, i.e. pre- 
dicts the emission in the central 100 AU diameter region. Since 
this area is completely inside the MIDI field-of-view, the model 
visibility can be directly compared to the observed visibilities. 
For the fits to the SED, however, we have to take into account 
that the very large beams of the far-infrared observations must 
clearly include emission from the surrounding large-scale cloud, 
on spatial scales of ~ 20000 AU (~ 0.1 pc) and larger, which 
is far outside the r = 50 AU model area. This large-scale cloud 
emission should not be confused with a possible circumstellar 
envelope, but represents the molecular clump and the surround- 
ing cloud in which IRS 1 i s embedded (see discu ssion below). 
The mm-maps presented bv lSchrever et alJ (Il997h show that the 
size of this cloud is about 1' (corresponding to ~ 55 000 AU or 
~0.27pc). 

The effects of this large-scale cloud are as follows: our MIDI 
measurements are insensitive to this large-scale emission, be- 
cause it is far more extended than the angular resolution (and 
also the field-of-view) of MIDI. Since the large-scale cloud 
emission is completely over-resolved for MIDI, it will not af- 
fect the observed visibilities. The observed fluxes, however, are 
expected to be strongly affected by the large-scale cloud. First, 
the cloud material in front of the embedded object IRS 1 will 
cause considerable extinction. This effect is represented by the 
foreground extinction we use in our modeling. The second ef- 
fect is that the cool dust in the large-scale cloud will produce 
far-infrared and mm-emission. Owing to the large beam sizes 
of the far-infrared and mm-observations, a large fraction of the 
observed fluxes in these beams will result from the large-scale 
cloud, and not from the central embedded object IRS 1. This ex- 
plains why our RADMC model fluxes at wavelengths ^ 20 //m 
are considerably lower than the observed fluxes. To approxi- 
mately include this large-scale cloud emission component in our 
model, we added to our SED model in Fig. |6] two additional 
blackbody components. For the temperature of the first com- 
ponen t, we used T = 55 K, as determined by ISchrever et alj 
(1199 7^ for the clump in which IRS 1 is embedded. For the sec- 
ond component, we used a temperature of 120 K, which is typi- 
cal of warm dust around ma ssive YSOs in so-called "hot-cores" 
dHerbst & van Dishoec3l2009.) . Fig |6] shows that with the addi- 



tion of these two components good agreement between the ob- 
served fluxes and the model can be achieved. 

An interesting question is whether this 120 K component 
could be considered as an envelope around the disk, as for one 
interpretation of hot cores as an infalling envelope undergoing 
an intense accretion phase (lOsorio et al.l 12009). The RADMC 
model of a disk without surrounding envelope provides a very 
good fit of both visibilities. Adding a spherical circumstellar 
envelope with a radius of < 250 AU would change the model 
visibilities significantly (by > 5%) and lead to disagreement be- 
tween the model and the data. On the other hand, a (more or less) 
homogeneous envelope with a radius of > 250 AU would be 
over-resolved for MIDI and not affect the visibilities. However, 
to reproduce the far-infrared fluxes in the SED, the radius of the 
120 K emitting region has to be at least > 400 AU (in the lim- 
iting case of optically thick blackbody emission). Owing to the 
limited MIDI field of view, the difference between the ISO and 
the MIDI spectrum could result from such an extended structure. 
The question of whether an extended circumstellar envelope is 
present can thus not be answered by the available data. Whether 
the material in such a hypothetical envelope would be gravita- 
tionally bound and ultimately accrete onto the circumstellar disk, 
remains another open question. 

4. Discussion 

We now want to discuss our findings presented in Sect. [3] We 
first compare the size of the mid-IR emission region as calcu- 
lated in Sect. 13. II to the corresponding sizes of similar objects. 
We then discuss the radiative transfer models. 



4.1. Comparison to MIR sizes of other massive YSOs 

The near-IR sizes of the disks of most YSOs (in particular 
Herbig Ae stars) sho w a tight correlation with the luminosity 
dMonnier et alJ l2005h because the emission in the NIR comes 
from regions close to the dust sublimation radius, which is sen- 
sitive to radiation from the central star and scales approximately 
with the square root of the luminosity. The mid-infrared emis- 
sion, in contrast, originates from a wider range of regions: the 
hot inner wall located at the dust sublimation radius, the surface 
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Fig. 6. Left: Comparison of the best-fit RADMC model (gray 40 m baseline, black 90 m baseline) for the visibilities of 
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a foreground extinction of 27 mag, and the black line the RADMC model with extinction and two additional blackbody components 
with temperatures of 55 K (yellow) and 120 K (pink). 



of an irradiated d i sk, an d the circumstellar envelope. As shown 
in iMonnier et al.l (12009 ). the size of the region emitting in the 
MIR therefore spans a much larger region for a given luminos- 
ity. For circumstellar disks, factors such as the flaring of the disk 
and the shadowing caused by the inner rim play an important 
role. Here, we wish to compare the MIR size of NGC 2264 IRS 1 
as found for the uniform disk model with the MIR sizes of other 
MYSOs, as well as with those of Herbig Be stars. 

Only a few other MYSOs have been studied with high 
angular resolution in the MIR until now. Amongst them are 
AFGL 2591, AFGL 2 136, AFGL 490 , and S 140 IRS 1 from the 
study of MIR sizes of Mon nier et al.l (2009), as we ll as W 33 A 
(Ide Wit et al. 2007. 2010) . and NGC 3603-IRS 9A (Vehoff et al.l 
201(31). The sample of Herbig Be stars was also taken from 
Monnier et al.l (12009*). The disk size vs. luminosity diagram of 
these objects and NGC 2264 IRS 1, as well as the size vs. stel- 
lar mass diagram is shown in Fig. |7] Parameters such as stellar 
mass, distance, and the numbers for disk size and luminosity can 
be found in Table lA.fl Both luminosity and mass of the central 
object seem to be correlated with the size of the MIR emission 
region. Owing to the factors considered above, objects with lu- 
minosities ranging between 1000 and 1 0,000 L© display a large 
scatter in their sizes (a factor of around 10). The same is true for 
the mass range of ~ 8 - 20 Mq. The size of the MIR emitting re- 
gion around NGC 2264 IRS 1 seems to be typical for their mass 
and luminosity. 



4.2. Discussion of the radiative transfer modeis 

Models of a circumstellar envelope only did not provide good 
fits, as they were unable to reproduce visibilities and NIR flux 
simultaneously. Furthermore, the estimated sizes of the dust dis- 
tribution (Sect. 13.1b are different for different position angles, 
which implies that there is an asymmetric dust distribution. 

The radiative transfer modeling performed with RADMC 
suggests that IRS 1 has an optical thick, but geometrically 
thin, flaring circumstellar disk with a mass of 0.1 Mq. This is 
supported by the results of the SED fitting process using the 
Robitaille grid. Not only the best-fit Robitaille disk model but 



also the ten best-fit models have parameters describing the cen- 
tral object and the circumstellar disk that compare very well 
with those of the best-fit model found using RADMC. However, 
the inclinations of the Robitaille models vary between 20° and 
70°. The temperature and mass of the central object are slightly 
higher in all of these disk models than the values used for the 
RADMC model. 

The two blackbody components producing the far-IR flux 
for our best-fit RADMC model could originate from a large sur- 
rounding cloud or envelope, respectively, that would with the 
large angular resolution of MIDI only be seen as foreground ex- 
tinction. This would agree with the large foreground extinction 
(> 20 mag) as used in the best-fit RADMC and derived from the 
ten best-fit Robitaille disk models. Furthermore, previous publi- 
cations have suggested that there is a large amount of foreground 
extinction with values between 20 mag (T hompson et al. 19981) 
and 35 mag ( Thompson & Tokunagal[l978 ). The inclination of 
the best-fit RADMC model is 30°, but owing to the poor uv- 
coverage this value is poorly constrained. The position angle of 
the jet-like feature seen in the near-IR is 20° (see Fig. [8]), and 
the position angle in our best-fit disk model is 40°. This does 
not match perfectly, but the value for the position angle is also 
poorly constrained. However, an overall geometric picture of a 
moderately inclined disk seems reasonable. 



5. Summary & conclusions 

We have presented MIDI mid-infrared interferometric observa- 
tions of the massive young stellar object NGC 2264 IRS 1. The 
observed visibilities provide no hint of multiplicity of the central 
source in the separation range from ~ 30 AU (resolution limit) to 
~ 230 AU (slit width limit). Geometric models (Sect. lTTl i sug- 
gest a size of ~ 30 - 60 AU for the mid-infrared emission region. 
This value is in good agreement with the mid-infrared sizes of 
similar objects. 

Our modeling shows that a simple temperature-gradient 
model for a circumstellar disk can reproduce the SED at mid- 
and far-IR wavelengths, but fails to reproduce the SED in the 
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ence of a large (r > 400 AU) spherical envelope around the 
disk cannot be excluded by the data. The mass of the circum- 
stellar disk is about 0.1 Mq. The data suggest an inclination an- 
gle of ~ 30° and a position angle of ~ 40° for the orientation 
of the disk. These values are consistent with an overall geomet- 
rical model based on the jet-like feature seen to the north-east 
of IRS 1. Compai-ing NGC 2264 IRS 1 to other YSOs, the size 
of its MIR emitting region seems to be typical of its luminos- 
ity and mass. However, the large scatter of sizes in this range of 
luminosities and masses points towards a wide variety of (disk) 
morphologies among these objects. 

More observational research, such as future MIR interfero- 
metric observations, will provide tighter constraints on the cir- 
cumstellar material around this interesting source and therefore 
will help us to clarify our understanding of the disks around 
high-mass YSOs. 
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Fig. 8. K-band image of NGC 2264 IRS 1 from lSchreveret alj 
(Il997h . We overlay an ellipse (white) that shows a disk seen un- 
der an inclination of 30° and a position angle of 40°. The size 
of the disk is not drawn to scale. The solid white line shows the 
orientation of the 40 m baseline, the dashed one the orientation 
of the 90 m baseline. 



near-IR and the observed visibilities; this model can therefore 
be ruled out. 

We have performed a detailed radiative transfer modeling 
of the observed SED and visibilities with the RADMC code. 
It suggests a scenario of a geometrically flat but optically 
thick circumste l lar du st disk. We also used the fitting tool by 
iRobitaille et alj (l2007h . which leads to similar disk models and 
therefore confirms the results obtained with RADMC, whereas 
envelope models found within the Robitaille grid did not pro- 
vide similarly good fits to SED and visibilities simultaneously. 
Although there is no observational evidence, the possible pres- 
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Table A.l. MIR-sizes of MYSOs and Herbig Be stai's 



Object 


Diameter 


Distance 


Luminosity 


Mass 


SpT 


Reference 




[AU] 


[pc] 


[Le] 


[Mo] 






NGC 3603 IRS 9A 


2177 


7000 


23-10" 


40 




1 


W33A 


115-230 


3800 


MO' 


10 




2 


V376 Cas 


23 


630 


43 




B5e 


3,4 


MWC 147 


19 


800 


1550 


6.6 


B6 


3,5,6 


RMon 


75 


800 


1400 


10.4 


BO 


3,7,8 


ZCMa 


61 


1000 


2400 


16 


B 


3, 9, 10 


MWC 297 


14 


250 


10600 


10 


B1.5 


3, 11, 12 


V1685 Cyg 


28 


980 


5890 




B3 


3,4 


MWC 361 


20 


440 


6600 




B2 


3,4 


LkHff 234 


27 


1000 


500 


5 


B7 


3, 13, 14 


AFGL 490 


100 


1000 


2-10^ 


8-10 


B2 


3, 15 


LkHo- 101 


49 


700 


2-6-10* 


15 


Be 


3, 16 


S 140 IRS 1 


122 


910 


2-10'' 






3, 17, 18 


AFGL 2136 


240 


2000 


2-10'* 


>10 




3, 19 


AFGL 2591 


118 


1000 


2-10'* 


10-15 




3,20 


M 17 SW IRS 1 


69-84 


2100 


S-IO^ 




BO 


21 


M8E-IR 


30-57 


1500 


2-10* 




BO 


22 


NGC 2264 IRS 1 


30-60 


913 


1.5-4.7-10^^ 


9.5 


B0-B2 


seeChaDt.miSH 



Refer ence s. (1) IVehoff e t al. (2010); ( 2) Ide Wit et al. (2007); (3) Mon nier et 
(12008'); (6) 'Hernandez et al. (2004); (7)'Murakawa et al. (2008); (8) MurakawJ 
(11) Ackeetal. (2008); (12) Drew et al. (1997); (13) Shevchenko & Yakubq; 
(16) Herbig et al. (2004); (17) Crampton & Fi sheil ( 119741) ; (18) ILester et al.l ' 
(21) ,FollertetaL(,201Q) ; (22) ,Linz et al. ,(2009.). 
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